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ABSTRACT 

We investigate properties of the reconnecting current layer in relativistic pair 
plasma reconnection. We found that the current layer self-regulates its thickness 
when the current layer runs out current carriers, and so relativistic reconnection 
retains a fast reconnection rate. Constructing a steady state Sweet-Parker model, 
we discuss conditions for the current sheet expansion. Based on the energy ar- 
gument, we conclude that the incompressible assumption is invalid in relativistic 
Sweet-Parker reconnection. The guide field cases are more incompressible than 
the anti-parallel cases, and we find a more significant current sheet expansion. 

Subject headings: magnetic fields — plasmas — relativity 



INTRODUCTION 



Magnetic reconnection (Vasyliunas 1975; Biskamp 2000 Birn & Priest 2007) is the 



driver of explosive events in space plasmas and laboratory experiments. Although it has 
extensively been studied ever since the 1950s, its detailed mechanism still remains unclear. 
The complexity of the reconnection problem is that small-scale physics in and around the 
reconnecting X-type region (or the diffusion region) can drastically change the system's 
global evolution. 

Magnetic reconnection has drawn attention in various high-energy astrophysical ap- 



plications (Coroniti 1990 Drenkhahn & Spruit 2002 Lyutikov 2003), where a relativistic 
extension of magnetic reconnection in an electron-positron pair plasma (e^) is considered. 
However, the physical mechanism of relativistic magnetic reconnection, as well that of the 
well-studied nonrelativistic reconnection, is far from being understood. For example, only a 
few researchers have developed fundamental models of relativistic steady state reconnection. 



Blackman & Field ( 1994 ) predicted that relativistic magnetic reconnection facilitates faster 
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energy conversion than nonrelativistic reconnection, due to the Lorentz contraction in the 



outflow region. Lyutikov & Uzdensky (2003) followed them, but Lyubarsky (2005) claimed 
that reconnection inflow cannot be so fast. He further worked on the generalization of rel- 
ativistic Petschek reconnection and argued that the guide field (out-of-plane field) opposes 
energy conversion into the plasma energy. 

On the other hand, self-consistent simulations have revealed important aspects of rel- 



ativistic magnetic reconnection. Several authors ( 


Zenitani & Hoshino 


2001 


2008 


Jaroschek et al. 


2004 


Bessho & Bhattacharjee 


2007 


) carried out particle- 



simulations of relativistic pair plasma reconnection, and they found that powerful particle 
acceleration occurs inside/near the diffusion region. The main accelerator is the reconnection 
electric field, which becomes strong enough to sustain relativistic outflow. When the system 
contains a guide field, relativistic magnetic reconnection also involves particle acceleration 



by the reconnection electric field (Zenitani & Hoshino 2005, 2008 Karlicky 2008). It seems 



that particle acceleration by the reconnection electric field is a common feature of relativis- 
tic magnetic reconnection. Furthermore, carrying out a relativistic magnetohydrodynamic 



(MHD) simulation, Watanabe & Yokoyama (2006) presented a Petschek structure in mildly- 
relativistic regime. Note that they assumed spatially limited resistivity, representing the 
physics inside the diffusion region. 

In an ion-electron plasma, it is known that a thin current layer appears near the X- 



point during magnetic reconnection (Horiuchi & Sato 1994). The main current carriers are 



electrons, and the thin electron current layer is a key region to understand the reconnection 
physics. So, significant efforts have been paid to the evolution of the electron diffusion region 



(Hesse & Winske 1998 Hesse et al. 1999). The general consequence will be that the current 



layer becomes thin as reconnection develops, typically into the electron meandering width 



(Horiuchi & Sato 


1994 


). On the other hand, there is an 


the current layer 


IFujimoto 


2006 


Daughton et al. 


2006 



may regulate the reconnection rate by changing the aspect ratio of the "diffusion" region. 

In a relativistic pair plasma, little attention has been paid to the detailed structure of the 
diffusion region, although it is a key region to drive magnetic reconnection and to accelerate 



high-energy particles. Recently, Hesse & Zenitani (2007) investigated the composition of 



the reconnection electric field, and they showed that the off-diagonal parts of the pressure 



tensor are important, similar to the nonrelativistic cases (Hesse et al. 1999, 2004 Bessho 



& Bhattacharjee 2007). However, a physical interpretation of the off-diagonal terms has 
not yet been established as well as the terms in a true decomposition of the relativistic 
pressure tensor. In the present paper we investigate the properties of the diffusion region in 
relativistic pair plasma reconnection from another viewpoint: the temporal development of 



-3- 



the reconnecting current layer. The unexpected new result is that the out-of-plane current 
layer adjusts itself by expanding its thickness. When the plasma drift velocity reaches an 
upper limit, typically on the order of the light speed, the current layer runs out current 
carriers. Then, the reconnection electric field becomes stronger to satisfy Ampere's law 
for the field reversal current. Consequently, the non-MHD diffusion region and the current 
layer become wider. Since it enhances the reconnection electric field and since it improves the 
aspect ratio of the dissipation region, relativistic reconnection retains a fast reconnection rate. 
A brief discussion of the steady state energy flow provides some insight into this process. We 
cannot employ the incompressible assumption in relativistic Sweet-Parker conditions any 
more. The guide field brings incompressibility, and therefore, we find a more significant 
current sheet expansion. 

The paper consists of the following sections. In §2 we describe our simulation setup 
and briefiy overview the system evolution. In §3 we extended our interpretation by using 
a steady-state Sweet-Parker model. In §4 we discuss the guide field extension. Finally, §5 
contains discussion and the summary. 

2. SIMULATION AND RESULTS 

We carry out two-dimensional PIC simulations in the x-z plane. As an initial cur- 
rent sheet configuration, we employ a relativistic extension of the Harris model. The 

magnetic field, plasma density, and plasma distribution functions are described by B = 
Botanh.{z/L)x + Boy, d{z) = (7/3^0) cosh.~'^ (z / L) and fs oc exp[— 7^{e — (3sUy}/T]. In 
the above equations, Bq is the magnitude of antiparallel magnetic field, Bq is the out-of- 
plane magnetic field, L is the typical thickness of the current sheet, no is the proper number 
density of plasmas in the current sheet, the subscript s denotes the species {p for positrons, 
e for electrons), (3p = —f3f. = (3 is the dimensionless drift velocity, 7^ is the Lorentz fac- 
tor for /3 (7^ — [\ — e is the particle energy, u is the relativistic four- velocity of 
u —\\. — {v/cf\~^/'^v, and T is the proper temperature including the Boltzmann constant. 
We set T — mc^ and (3 — 0.3, respectively. In addition, a uniform background plasma 
is added to the system. Its number density and temperature are 7165/(7/3^0) = 5% and 
Thg/mc^ = 0.1, respectively. This background temperature is cold enough that the back- 
ground plasma energy density is ~ 1.2niygmc^ per species, including the rest-mass energy. 
The plasma satisfies the pressure balance condition Bf^/Sn = 2noT in this equilibrium. The 
Debyc length is [T / (An'j'^noq^)]^^'^ = 0.3L, where q is the charge. The typical Larmor radius 
is c{qBo/'ymc)~^ ~ 0.2L. 

The system consists of 1536(a;) x 768(2;) cells, and the typical scale of the current sheet 
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L is set to 10 cells. We consider periodic boundaries in the x-direction, and the boundaries 
are located at x = ±76. 8L. We also consider periodic boundaries in the z-direction, but 
we divide the simulation domain into two subdomains. The first half (the main simulation 
domain) is located in —19.2 < z/L < 19.2, and the other half (the subsimulation domain) is 
located in 19.2 < z/L < 57.6. In the subdomain, we assume an oppositely directed Harris 
model (e.g. B = — i?otanh[z/L — 38.4]i), so that physical properties in the sub domain 
are smoothly connected to those in the main domain. We use 7.6 x 10^ superparticles in 
this simulation. One cell contains 6.4 x 10^ particles at the center of the current sheet. The 
detailed parameters are presented in Table [T} During the very early stage of the simulation, 
we set a small driving force to trigger magnetic reconnection near the center of the simulation 
domain. This trigger field smoothly varies in time, and it is applied until t ~ 12rc, where 
Tc = L/c is the light transit time, while we discuss the physics of reconnection in the late 
stage of SOtc ~ 90tc. The total energy is conserved within an error of 0.1% throughout the 
simulation run, after the initial trigger force vanishes. In addition, we carry out four runs 
with different parameters of n^g / {'y prio) and Bq/Bq. Their parameters are shown in Table [ij 
These configurations are similar to previous simulations (Zenitani & Hoshino 2007, 2008[ ). 



Here we describe the evolution of the first run (run 1). After the initial impact dis- 
appears, magnetic field lines start to reconnect around the center of the main simula- 
tion domain. The panels of Figure [T] present the snapshots of the reconnecting region at 
t/Tc = 70. The X-point is located at the left border {x = 0), while the simulation domain 
is —76.8 < x/L < 76.8. The top panel shows magnetic field lines and the y-current (out-of- 
plane component of the current) density. The bottom panel shows magnetic field lines and 
the ?/-component of the plasma average velocity. Notice that the X-point is located around 
X ~ 0. At this stage, magnetic reconnection is well developed: the plasma outflow velocity 
{{vx)) is up to ~ 0.7c around the outflow region. 

Panels in Figure [2] present the plasma number density 2n, the current density Jy, average 
drift velocity (vy), and reconnection electric field Ey along the inflow line x = 0. These values 
are averaged over —1/4 < x/L < 1/4 in order to reduce noise. The y-current is carried by 
the drift motion of positrons (+|/-direction) and electrons (— y-direction) inside the current 
sheet, Jy ~ 2d{z){vy). We present three characteristic stages at t/Tc = 30, 50, and 70. Shortly 
before t/Tc = 30, reconnection is about to begin: {vy) and Ey start to increase around the 
X-type region, while plasma density starts to decrease there, because plasmas flow away 
into the reconnection outflow region. At t/Tf. = 50, the reconnecting current layer looks 
"thin" around the reconnecting region, and then the plasma inflow and outflow structures 
are developed. The {vy) goes up to ~ 0.7c around the center. 

After t/Tc = 50, the current layer shows an unexpected evolution which is not explained 
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Fig. 1. — Magnetic field lines and y-component of (a) the current density and [b) the averai 
velocity at t/xc = 70 in run 1. 
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Fig. 2. — Physical properties along the inflow region at x/L = in run 1: (a) normalized 
plasma density, {h) y-component of the electric current, (c) y-component of the plasma 
average velocity, and {d) reconnection electric field Ey. The red lines show the relevant 
{—V X B)y electric fields. 
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by present reconnection theories. Previous studies show current layers, which are substan- 
tially thinner than the initial current sheet. However, in this case, the current sheet starts to 
expand. The change may look small in Figure [H but we can clearly recognize the expansion 
signature in Figures |2|) and[2]c. The plasma velocity {vy) and the current density Jy increase 
outside the initial current sheet. In fact, Jy outside the current sheet is stronger than in 
the initial state, which was almost zero. Figure [sja {black line) presents the time evolution 
of the current sheet thickness measured by the scale height or the relative amplitude of the 
current density \Jy\- During the preflight stage, the current sheet becomes slightly thinner 
than the initial Harris sheet thickness, but it starts to expand, and finally it becomes ~ 1.9 
times wider after I/tc = 70. 

Figure |3^ shows time evolution of the reconnection rate, represented by the electric 
field Ey at the X-point. In this two-dimensional configuration, the reconnection rate Ey/Bo 
immediately indicates the efficiency of magnetic energy conversion. The rate is often nor- 
malized by the reconnection outfiow speed in nonrelativistic studies, but we normalize the 
rate by c, since the outfiow is on the order of ~ 0.7c. So, although this rate may be slightly 
underestimated, we obtain a relatively fast reconnection rate of 0.15-0.2. This rate is faster 
than the well-known nonrelativistic reconnection rates of ~ 0.1, due to the E'y-enhancement 
along with the current sheet expansion. 

Next, let us interpret the physics of the current sheet expansion. In this simulation, we 
use a large density ratio of Harris sheet plasmas to the background plasmas, n^g/ (7/3^0) = 5%. 
Once reconnection starts, dense Harris sheet plasmas escape into the outfiow region, while 
only less dense background plasma enters the reconnecting region. Therefore, the plasma 
density continues to decrease as seen in Figure |2]a. The density finally decreases to about 
the background density. Consequently, the system starts to run out the y-current Jy for the 
field reversal around the diffusion region. Recall the restriction of Ampere's law: the system 
needs the out-of-plane current Jy, 

dz c c dt 

By some mechanism, the system needs to maintain the field reversal. 

The first workaround is to enhance the plasma current Jy = 2dcs{vy), where dcs is the 
number density in the current sheet. In the nonrelativistic regime, the plasma current Jy is 
simply enhanced by particle acceleration into the ±?/-directions. The displacement current 
works indirectly, because the enhanced electric field Ey accelerates particles. In fact, at 
t/Tc < 50, the reconnection electric field Ey and the plasma y-velocity (vy) are significant 
only around the center {z ~ 0; Figs. [2]c and[2]fi). However, in a relativistic plasma, the 
plasma current has an upper limit of \J\ < 2dcsC. When the velocity becomes on the order 
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Fig. 3. — (a) Temporal evolution of the current sheet thickness at x/L = 0, measured by 
the scale height of the current density (log J^). {b) Temporal evolution of the reconnection 
rate (Ey/Bo) in run 1. 
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of c, plasmas cannot carry any more current. 

The next workaround is to directly use the displacement current {l/c){dE/dt). Indeed, 
in this case, the plasma current explains only 80% of the required current and the displace- 
ment current accounts for the other 20% around the turning phase of I/tc ~ 50. However, 
the system cannot use the displacement current to maintain a steady or quasi-steady struc- 
ture. This workaround is temporal, and then the system goes into the other phase around 
t/n ~ 50. 

Because of the effects of the displacement current term, the current layer starts to 
expand by the derivative of the generated electric field. Because of the displacement current, 
the electric field Ey becomes even stronger across the current layer. This reduces magnetic 
fields outside the current layer. 

Note that -B^ > and {dEy/dz) < on the upper side {z > 0), while -B^ < and [dEy/dz) > 
on the lower side {z < 0). Consequently, the MHD frozen-in condition breaks down over the 
wider spatial region around the neutral sheet. The red lines in Figure |2](i present (—vxB)y 
at the relevant stages. We recognize that the non-MHD diffusion region of Ey ^ {—v x B)y 
becomes wider during t/Tc = 50 ^ 70. Obviously, the non-MHD region corresponds to the 
plasma drift region of {vy) > 0. Since plasmas are free from the frozen-in restriction within 
the diffusion region, they are accelerated into the ±y directions by the enhanced electric field 
Ey, and then they start to carry the ?/-current. At t/r^ = 70 we can recognize that the outer 
regions around z/L ±2 carry the y-current (Figs. |2]6 and|2]c) and become a part of the 
diffusion region (Fig. ^fl). The enlarged diffusion region involves more plasmas outside the 
initial current sheet, and then they are responsible for the required Jy for the field reversal. 

In Figure |3] we can compare the trend of the current sheet thickness and the reconnection 
electric field. The reconnection electric rate Ey increases first and then it stops before 
t/Tc ~ 60. The current sheet thickness slightly delays to the reconnection electric field, and 
it starts to expand. The current sheet remains nearly constant after t/Tc = 70-80. 



3. ANALYTIC THEORY 

In order to understand the current sheet expansion, we consider a steady state Sweet- 
Parker current sheet model as presented in Figure |4| This model is too simple to discuss the 
detailed properties of simulation runs, but sufficient to understand the physics. In addition, 
as long as we consider a mildly relativistic regime of our scope, the current layer has a simple 
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planar structure around the entire reconnecting region (Fig. [T]). In Figure |4] the subscripts 
"in" and "out" denote the physical properties in the inflow region and in the outflow region, 
respectively. The current sheet width and height are 2Lm and 2Lout- In the steady state, 
conservation of magnetic flux, the number density continuity and energy budget in the model 
current sheet are written as 

BinVin = BoutVout (3) 
'^in^inLin'Vin '^out^outLout'^out (4) 

/ B^ \ / \ 

+ LinVin = [2-fl^^Wout + LoutVout, (5) 

where w = e + p is the plasma enthalpy density, the sum of the plasma internal energy (e) 
and the proper pressure (p). We introduce the factor of 2 in equation |5] so that we can 
consider both electrons and positrons. In the x-z plane, we can assume that the two species 
move together. By using the ratio of Poynting flux to the particle energy flux 

B' 

471(272^;)' ^ ' 

equation [5] can also be written as 

Note that the relativistic Alfven velocity is Va = [o"/(l + cr)Y^'^c and the relevant Lorentz 
factor is 7a = (1 + cr)^/^^ ^ cold plasma limit, ain is roughly reciprocal to the background 
plasma density (n.bg/'ypno). Considering the initial equilibrium condition B^/8tt = 2noT and 
the inflow plasma number density 'jinnin ~ Ubg, some algebra yields 

2noT 2 7/3^0 T 



0"r 



In the case of nonrelativistic antiparallel reconnection, in the limit of {vm/vout) = <C 1 
these conditions lead to 

dinLinVin doutLout'Vout (9) 

R2 2 

(^)L,„t;,„ ~ {2douJ^)LoutVout. (10) 

Note that we use the simulation frame density d = 'yn here. We obtain a familiar estimate 
of 



Bf 

Vout ~ A / 7 , 
ATrma. 



--V2 Va^u. ill] 
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Fig. 4. — Simple model for a Sweet-Parker current sheet. Inflow properties and outflow 
properties are represented by subscripts "in" and "out," respectively. 

Notice that the Alfven speed in an electron positron plasma is Va = B /[4:Tc{2md)Y^'^ in the 
nonrelativistic limit. 

Let us consider the relativistic case. For simplicity, we assume that incoming flow is 
flux dominated: the magnetic energy carries most of the incoming energy. So, we drop 
Win from the left-hand side of equation [5j This approximation is appropriate for our low- 
background density runs, because we set low density background plasmas {nhg/[ji3nQ] -C 1) 
in a nonrelativistic temperature {Tbg <^ mc^). In addition, we approximate the outflow 
enthalpy term in the relativistic pressure limit. 

Wont = {eout + Pout) ~ (uoutmc^ + ~ 
where F ~ 4/3 is the adiabatic index. In the outflow region, the pressure balance across the 

2pout (13) 



current sheet (Lyubarsky 2005) yields 

Stt 



Using equations and 13 and introducing u = Vin/vout, we obtain 



(^-^^-^)4:PoutLinVin ~ {Sj'^^^ + 4:U^)poutLoutVout (14) 
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The number density yields 



'yout'^out 



Lout'^out 



1 + CTi 



{Hut + 



In the hmit of -C 1 and 3> 1, we obtain 

yout^out 



27, 



2 

out 



(15) 



(16) 



in the simulation frame density, or {uout/nin) ~ '^yinlout in the proper density. This analysis 
tells us that the relativistic steady state reconnection is not incompressible in both the simu- 
lation frame and the proper frames. Previous reconnection models ( Blackman &: FieTd||1994 



Lyutikov 2003) employ an incompressible assumption in the proper frames. Our analysis 



shows that the relativistic reconnection model should be re-constructed by taking compress- 
ibility into account. Incompressible models are helpful when and only when '-y^ut ~ 1- 



We can confirm the above argument by using Lyubarsky (]2005)'s analysis. For two 

(17) 



fluids, the x-momentum conservation along the outflow line yields 

1 



O ^ ( llut'U^outvlut . 
2—1 5 VPaut 



dx 



~ -JyBout 
C 



Substituting (d/dx) ~ 1/i^m and Jy ~ cB in /{Arc Lout) and using equation 12, we obtain 



Hut^, 



2 

out '^out 



+ 1 ]Pout 



Bi 



2 ^T^Lout 



-B. 



out 



Combining equations |4, 13, and 18, we obtain 

Hut{Vout/cf + 1 



yout'^out 

yin^in 



Lout^out 



Hut 



(19) 



Next, let us consider the current sheet thickness. In the steady state, the current in the 
central current sheet satisfies 

cB 



2qdcsn 



out 



(20) 



where dcs is the typical number density in the current sheet. Also, the initial Harris condition 
satisfies 



cBq cBir, 



A'kL AttL 
The ratio of the current sheet thickness yields 

^out 



2g(7/3^o)/3c. 



(21) 



(22) 
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We can evaluate the right-hand side by considering the inflow number density ~ linnm, a 
central current sheet density 'jinnin < dcs < loutnout, and the compressible factor of equation 



16 



(3c (7/3^0) 1 13c 

^ rs_/ _ 

ribg dcsVy ~ nbg Joutnout Vy Ubg 27^„t Vy 



Then, equation 22 can be rewritten as 

Lout > (7/3^0) P ( c 



(23) 



(24) 



We recognize that the following condition is sufficient for the current sheet expansion, 

where Z) is a discriminant term. Note that D also indicates the minimum relative thickness 
of the steady state current sheet. Lout ^ DL. 

Since we do not have a perfect theory on 7out, we refer to the simulation results. In 
the case of run 1 (t/xc = 70), the plasma x- velocity Vout and the electromagnetic field 
do not satisfy the frozen condition inside x < lOL, and so we regard the region as the 
diffusion region. Inside the diffusion region, the Lorentz factor of the averaged flow 7^ = 
[1 — ((fi-)^ + (vy)"^) / c^]~^ ~ 2 stays constant in the diffusion region. The density stays small 
dcs ^ 2n6g and it gradually goes up in the outflow direction. Taking the results into account 



(e.g. 27^^^ ~ 4), we flnd that the condition (eq. 25) is already satisfled, D = 1.5 > 1. The 
observed thickness ~ 1.9L is larger than 1.5L. For comparison, we carried out the other 
run with dense background plasmas (run 2; Ubg/'j/snQ = 20%), too. In this case, the current 
sheet expansion is less signiflcant (thin line in Fig. [3]a). Although the observed thickness 
~ 1.5L is rather larger than D ~ 0.6 primary due to the slow y-velocity of (vy) ~ 0.5c, it 
is reasonable that the current sheet expansion is less signiflcant, because the inflow delivers 
more current carriers than in run 1. 



4. GUIDE FIELD CASE 

We also study the current sheet problem in the guide fleld case. Snapshots of run 
3 with strong guide fleld Bg/Bq = 1.5 are presented in Figure |5} They are at the well- 
developed stage of t/Tc = 140. Note that the guide fleld reconnection evolves slower than 
the anti-paralell reconnection. The typical outflow speed is Vout ~ 0.3c. The structure of 
the reconnection region looks more complicated than the anti-parallel case. For example. 
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we can see inclined current layers along the separatrix (Fig. Isla). There are also in-plane 



electric fields (E^^Ez) due to the charge separation of positrons and electrons (Karlicky 



2008 Zenitani & Hoshino 2008). However, in order to discuss a zeroth-order structure of 
the diffusion region, a similar Sweet-Parker approximation would be plausible. In fact, the 
^/-velocity structure, which is a good indicator of the diffusion region structure, looks like 
a rectangular structure (Fig. |5]6). Panels in Figure [6] present the plasma properties along 
the inflow line x = 0, similar to those of Figure |2} Surprisingly, the plasma number density 
in the current sheet is slightly lower than or almost the same as the inflow density at the 
well-developed stage {t/Tc = 140; Fig. [6]a). Plasma seems to be incompressible unlike the 
anti-parallel case. The plasma average y- velocity is even faster: (vy) ~ 0.9c at t/rc = 80 and 
{vy) ~ 0.97c at t/Tc = 140 at the center. This is probably because the guide field confines 
particle motion into the y-direction. Importantly, the current sheet looks much thicker than 
the anti-parallel case (Figs. |5]a, |6]6, and[6]c). 

Next, let us consider a similar Sweet-Parker model. We set the guide field in the infiow 
region to -Bcm- The energy fiux and pressure balance condition are modified in the following 

way, 

r2 I r2 r2 I r2 



f^in^^G,in\T fc^ 2 , ^ out ^ ^ G ,out\ j /nr\ 
y jLinVin = {2'J^^^Wout ^ jLoutVout (26) 

+ 2pouX (27) 



r2 I r2 ^ r2 

^in 'T' ■^G,in ( -^Gfiut 



Stt V 87r 

Hereafter, we dropped the cx-related terms since [(1 + o"i„)/crj„] ~ 1. The guide field cases 
are more fiux dominated than the anti-parallel case due to the guide field magnetic energy 
(Table [T]). The effects of in-plane electric fields are included in the Poynting fiux terms, as 
guide field contributions. In addition, fiux conservation of the out-of-plane field yields 

Bg Bg (28) 



Combining equations [T2| |26} and [27j, we obtain 

LinVin ^llutVout + {,Bl^^ + 5| „„J / An 

_ 

^out'^out Apout + {Bl.)lAn 

This immediately tells us that 



(29) 



and then equation 28 yields 



7 Z i> I'^Uj 



Bgm ^ Bcout- (31) 



- 15 - 




10 20 30 40 50 60 




10 20 30 40 50 60 



Fig. 5. — Magnetic field lines and y-component of (a) the current density and (6) the average 
velocity at t/xc = 140 in run 3. 
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Fig. 6. — Physical properties along the inflow region aX, x/L — Q (a) normalized plasma 
density, (&) y-component of the electric current, and (c) ^/-component of the plasma average 
velocity 
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So, the guide field is compressed in the outfiow region. For guide field larger than the anti- 
parallel component Bin ^ Bcm, more energy is stored in Bcout in the outfiow region. As 
the guide field becomes even stronger, the guide field term becomes the largest contributor 



in equation 29 and then equation 30 becomes closer to the unity. This means that the guide 
field brings incompressibility in the simulation frame. In the strong guide field limit, we 
employ an incompressible condition 

dcs ~ loutnout (32) 



and then we obtain a condition similar to equation 25 for the strong guide field limit 



Do = > 1, (33) 

Ubg 

where Dq is the discriminant term or the minimum relative thickness. It is significant that 
Dg only depends on the initial configuration parameters, while the antiparallel counterpart 



D (eq. 25) contains a variable 'jout, which comes from the plasma compressibility effect (eq. 



22). As a result, the sufficient condition is more easily satisfied than in the anti-parallel case. 



so that the current sheet expansion is more significant. 

The temporal evolution of the current layer thickness in the guide field runs is presented 
by the dashed line in Figure [3]a. In the guide field cases, time is re-arranged by some offsets 
(At = GOtc) due to the late onset of reconnection, so that we can directly compare their 
thickness in the expanding phase. In the guide field cases, the current sheet originally 
becomes thinner during the long preonset stage of reconnection. In the case of run 3, the 
outfiow guide field is Bcout ~ 1-7-Bo and so the system seems to be sufficiently incompressible. 



The sufficient condition (eq. 33) yields Dq = 6 > 1, which predicts that the current 
sheet expands up to a factor of 6 in the strong guide field limit. Considering the weak 
compressibility, this is consistent with the observed thickness of 4L-5L (Fig. [3]). In the 
case of run 4, the current sheet expands slowly (dotted line in Fig. |3]a), due to the slowest 
evolution of reconnection. The condition is Dq = 1.5 > 1 in the guide field limit, which is 
roughly consistent with the asymptotic thickness of the current sheet, ~ 1.3L. The parameter 
Dq may be a good approximation of the thickness in the guide field cases. One reason is 
that plasmas are nearly incompressible. The other reason is that the plasma y-velocity is 
closer to the light speed {vy) ~ 0.9c, because their motions are threaded by the guide field. 



5. DISCUSSION AND SUMMARY 



We demonstrated that the relativistic reconnecting current sheet self-regulates its thick- 
ness, so that it facilitates the fast reconnection rate. In order to study the aspect ratio of the 
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diffusion region, much attention has been paid to the length of the diffusion region. However, 
it is found that the relativistic reconnecting layer changes the aspect ratio by expanding its 



thickness, and that we obtain a higher reconnection rate. Zenitani & Hoshino (2001) briefly 



reported the current sheet expansion and argued that the relativistic mass increases the me- 
andering width. This explains some aspects, but the situation is more complicated. Indeed, 
due to the enhanced reconnection electric field Ey, some incoming plasmas start (relativis- 
tic) meandering motion around the neutral plane (2; = 0), while others are directly driven 
to the ±y-directions without reaching the neutral plane. The current sheet expansion arises 



from the net effect that contains all particle motions there. Bessho & Bhattacharjee (2007) 



reported the faster reconnection rate in the lower density background condition. That trend 
can be consistently explained by our theory. Based on a relativistic Sweet-Parker analysis, 
we developed physical interpretations and derived some critical conditions. It is found that 
the current sheet is likely to expand when the reconnection inflow is flux dominated, and it 
further expands when the guide fleld is strong. Although the guide fleld effect is unclear in 
the dense background runs (runs 2 and 4) , low-density cases exhibited a signiflcant difference 
(runs 1 and 3). In other words, in the high-cr (low-density) limit of astrophysical interest, 
the guide fleld may play a crucial role in the reconnection topology by expanding the current 
sheet. 

We also found that the relativistic Sweet-Parker model should be treated compressibly. 
In fact, our simulation snapshot shows the density ratio of > 2 in the steady state (e.g. Fig. 
2b). Previous authors assumed plasma incompressibility (Blackman & Field||1994 Lyutikov 



& Uzdensky 2003), but these models should be re-constructed under the compressible con- 



dition, especially when the outflow Lorentz factor can be large {'^out S> 1), unless the guide 
fleld is strong enough. The relativistic Petschek models ( Blackman fc Field||1994 Lyubarsky 
2005| [Tolstykh et al.|[2007l ) are out of the scope of the present study, but our analysis on the 



guide fleld reconnection supports Lyubarsky (2005)'s argument that the guide fleld changes 



the dynamics of reconnection. In our case, the guide fleld introduces the incompressibility, 
and then the current sheet is more likely to self-regulate. 

We flnd another new result regarding the outflow speed. In the nonrelativistic regime, 
the reconnection outflow speed is often approximated by the inflow Alfven speed. If this 



can be applied to the relativistic cases 7^^^^ ~ l\in ~ (1 + o"m) ~ cTm, equation |8] implies 
that the critical condition (eq. 25) is insensitive to {"^ i^tiq) / Ubg or am- However, in our 
simulations, the obtained outflow Lorentz factor -^out is substantially slower than that of the 
inflow Alfven speed (7A,m)- The outflow velocity is insensitive or weakly sensitive to the 
inflow flux (Tin, and therefore, the current sheet expansion does depend on {'y/3no)/nbg. It is 
reasonable that the current sheet expansion is more signiflcant when reconnection inflow is 
more flux dominated, because it contains fewer current carriers. There are several reasons 
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why reconnection outflow is "slow." One reason is that reconnection is unsteady. The 
reconnection jet slows down because it pushes away the preexisting plasmas at the outflow 
region in unsteady simulations. This will be checked by the PIC simulations with open 
boundary conditions. Another is the relativistic pressure effect. In the relativistic regime, 



a relativistically high pressure increases the plasma enthalpy (eq. 12), which works as an 



effective inertia, and then bulk speed becomes slower. In order to maintain the current sheet 



against the strong magnetic pressure (eq. 13), a relativistically high pressure is required. We 
can see this from a calculation which neglects pressure effects. Without the pressure terms, 
the energy and number conservation (eqs. Ill and [5| read 



^1 + CTi^) „ (34) 
'^I'out^ out _ itut^out _ 'Jout LinVin 



From the equations, we would obtain an unrealistic super- Alfvenic outflow, -^out ~ (1 + 



(^in)lin 3> 7A,m (Lyutikov & Uzdensky 2003). As a result, we conclude that the pressure 
effect is critical. 

Compared with the one-fluid approximation, our model deals with both electron fluid 
and positron fluid independently. Therefore, the number/energy budget conditions are more 
reliable, and we can also discuss the ?/- velocity \vy\ < c. Instead, the difficulty is that several 
parameters {dcs, Jout, Vout and Vy) still have freedom in the current sheet or in the outflow 
region. The electron/positron fluids with high Lorentz factor 'jout may flow toward the x- 
direction or toward the ±?/-directions. The reasonable criteria {Vy + vl^^) < in the outflow 
region will improve the model. 



It is reasonable that the conditions (eqs. 25 and 33) are reciprocal to the drift speed 
parameter f5. When the initial drift speed is fast, the current sheet is relatively thin, so that it 
is more likely to expand. However, in the relativistic speed limit of /5 — » 1, the current sheet 
is likely to expand before the large-scale reconnection evolves due to the various instabilities 
(e.g. the tearing instability or the drift kink instability). In the vacuum inflow limit of 
nbg — ^ 0, eventually magnetic reconnection will disappear because there are no plasmas to 
carry the steady current. The induced electric fleld simply travels away as an electromagnetic 
wave, without accelerating plasmas. We do not know whether the transition is smooth 
or drastic, from the self-regulated current sheet to the quiet vacuum condition. The self- 
regulation of the current sheet will be found in relativistic reconnection in an ion-electron 
plasma, too. When the electron drift speed increases to the order of c, the electron current 
layer will similarly expand, while it is unclear that the enhanced electric fleld accelerates 
ions. When the broadened electron layer becomes comparable to the ion layer, ions and 
electrons more effectively interact with each other, and then it may affect the reconnection 
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structure once again. 

In summary, we demonstrated that the relativistic reconnecting current sheet self- 
regulates its thickness, in order to carry the required current. A simple Sweet-Parker analysis 
shows that the current sheet expansion is dominant in the flux-dominated inflow and that the 
guide field enhances the current sheet expansion. It is also noteworthy that the relativistic 
Sweet-Parker reconnection is compressible and that it is incompressible only in the limit of 
the strong guide field. 



Table 1: List of Simulation runs 



Run 


1 


2 


3 


4 


T/mc^ 


1 


1 


1 


1 










1.5 


1.5 




0.05 


0.2 


0.05 


0.2 




~ 36 


~ 9 


~ 115 


~ 30 


VoutI C 


~ 0.7 


~ 0.4 


~ 0.3 


~ 0.3 


i<mt 


- 2 


- 1.3 


~ 3-4 


~ 1.2 
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